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ABSTRACT 


In Paper I, we derived equations for the perturbations on an artificial satellite 
when the motion of the satellite is commensurable with that of the earth. Tliis was 
done by first selecting the tesseral harmonics that contribute the most to the per- 
turbations and then by applying Hori's metliod by use of Lie series. Here, we intro- 
duce some modifications to the perturbations, which now result in better agreement 
with numerical integration. 
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ON THE TESSERAL-HARMONICS RESONANCE 
PROBLEM IN ARTIFICIAL-SATELLITE THEORY’^ 

Part II 

Barbara A. Romanowicsc^ 


1. INTRODUCTION 

If tlie gravitational potential of the earth is expanded in terms of Legendre poly- 
nomials and functions, then in order to obtain a good approximation in tine determina- 
tion of the orbit of an artificial satellite, it is usually sufficient to consider the zonal, 
longitude-free terms of the expansion. However’, the influence of the tesseral terms 
becomes important when the moan motion of the satellite and the rate of rotation of 
the earth around i ts axiS are in a simple ratio; this is called tosseral-harmonics 
resonance. 

In Paper I (Romanovdez, 1975), a theory was derived to compute analytically the 
perturbations on the moi.ion of an artificial satellite in the case of tesseral-harmonics 
resonance. In a comparison with a numerical integration, however, some discrepan- 
cies appeared in the analytical method. After some modifications wore made to the 
perturbation expressions, which we derive herein, the theory has proved to agree much 
better with numerical integration. 


* 

This work was supported in part by Grant NGR 09-015-002 from the National Aero- 
nautics and Space Administration. 
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2. NEW EXPRESSIONS OP THE PERTURBATIONS 

2. 1 Perturbing Fimction and Equations of Motion 

The gravitational potential of the earth at exterior points can be expressed as 
follows [cf. equations (1) and (2) of Paper I]: 


Ji 


U(r,0,X)-^ 


‘-f ^n(f 


n-^2 


XI r- 2 m--l 


n 


m \ r / n, m 


P (cos 0) cos m(\-X, 

( 1 ) 


h,m 


where (r, 0, X) are spherical polar coordinates relative to the center of mass of the earth, 

the tuiis of rotation being the pole of coordinates; p is the gravitational constant G times 

til 

the mass of the earth; R is the mean equatorial radius of the earth; is the n 

Legendre polynomial; P (z) are associated Legendre functions; 

IX^ ill 


Pn,m(^') 


d?- 


and J and ^ are dimensionless coefficients I’clated to the normalized coefficients 
n n,m 

^n, m’ ^n, m 


‘^n,m'"' '^n,m \/2(2nM) (n-m)! /(n i m)! , m v*.- 0 , 


-Jn=^Cn^0= V^^n, 


( 2 ) 


■^n, m ^ ^ cos mX S„ sin mX . 


n, m' n, m 


n,m 


In choosing the perturbing function, we talce the leading zonal harmonic (i.e. , 
the one containing Jg) and we select those tesseral harmonics that, in the particular 
resonance being studied, give rise to long-period effects with the largest amplitude. 
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The general tesseral harmonic can bo developed in toimiB of the osculallng elliptic 
elements (a, o, I, il, w, M), referred to tlie equator of the earth, in the form (Kaula, 
19G6); 


nfl 


n, m (ll)(a ) ‘^n, m j n, ni, ^n, p, q^*^^ 


p-n 


I 


(n-2p)u) ^ (n~2pf q)M 


P“0 


/ 




( 3 ) 


where is the angaiar velocity of the rotation of the earth, t is the time, and F (I) 
find q(e) are, respectively, the inclination jmd eccentricity functions as defined 
by Kaula (196G). The functions G (e) are of order |q| in eccentricity. 


Resonance occurs when a pair (a, (3) of mutually prime integers exists such that 
the satellite performs p nodal periods while the earth rotates a times relative to die 
processing satellite's orbit plane. This can be eijqpressed by 


a(u +M) p(v- fi) , (4) 

where w, M, and are the rates of change with time of w, M, and S2, respectively. 
The corresponding slowly varying ai’guments are of the form 

- f const , 

1 a, p ’ 


where 


a(co + M) + P(«-v/t) 

is called the resonant variable and - 1, 2, 3, . . . . 

Considering that to is a small quantity and that the general argument in a tesseral 
harmonic is 
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4> - (n-2p)^ t (n* 2p'ti)M f _) , 

IX* 

we shall keep those tesseral harmonics that contain arguments such that 
n - 2p (' q * kja , 

ki • « ft ft 

m ‘ kjP , ‘ 


(r>) 


Since n ‘'m and since lower order tesseral harmonics are bound to have larger 
effects because of tlie factor (Il/r)^, whore R/r < 1, it can be assumed that it is suffi- 
cient to consider only tlie cases kj^ = 1, 2, and 3ft In a general manner, the tesseral 
hairmonlcs to be kept are 

\j,p+kQ,k^p ^ 

where « i. 2, 3, . , , and for each kp tlio index takes vulucs 0, i, 2, , . . . Wo 

then write each V in tlie form 

n, m 

V = V' + R , 
n, m n, m n, m ’ 

where R^ the "resonant" part of the tesseral harmonic, consists of all the terms 
that contain the resonant arguments, obtained by solving e(iuatims (0) for p and q 
[cfft equations (6), Paper I]. Thus, the resonant paxt of the tossox'al hamxonic V 

11 ^ 111 
is 


pkj-ko 




where 


' u 




( 6 ) 


q^^^ kj(a-p) - kg 2x 


and 

The residua’ part of is whose effects are much smaller than those from 

^^n, m* 

In Paper I, we then defined D(k^) and 0j^(k2) by 
where 


D(kp cos 20^(kj) ^ \ (j( k„,k,p 


the sum over kQ(kp meiming that we have taken into account all the values of k^ when k^^ 
has a given value. In order to obtain D(kj^) and 0|(k|), we can write 


S(ltpkQ) 




n,m 


n f 1 


a 


n=k^pJkQ , m-‘-kjP , 


A(k,,krt,x)"F (1) G 
' P 0’ ' n,m,x’ ' ] 


n,x,qjj 


(e) 


(7) 


X'n 

D(kj^) exp 2i0^(kj) - E ^ V -^(l^p^^O'X) exj) - qw - 

kpOc^) x-0 ’ 


If we write 

\,ko"^^"#\p4-ko,k^p 
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and 


2aj(kj) - kj(a-(3)w f i.'(kj) , 


then (see Appendix A) 


x--n 


D(kj) oxp l!/(k^) " E «(lvi, Hq) ^ A(k J, k^, X) t'xp i [(kg - 2x) w - Aj^ . (8) 

ko(kj) 0 


Finally, the po.ential in which the Batollite moves is 




m 


( 8 ) 


the sums being taken over ail the tessoral harmonic.'^ selected. 


If we consider the canonical set of Delatmay variables, 


./jia , yiia(l-e2) , cos 1 , 


"D 




*D D 

hjj " S2 , 


( 10 ) 


then the Hamiltonian of the problem is 


r,; -li 4 T<’ , 4- F 

ot2 1 2 ' 


(11) 


where 


1^2 

Fi= ’ 


2 ” ^ ^ ^n, m ‘ 
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To avoid tlio depondonce m time of the Hamiltonian, we perform a canonical trana- 
formatlon, bo that the new variablen Icf, equatlona (12) of Paper I] are 


L “ Lijj ) 
G - Gr 


'D 


II -• II 


D 


K-Kp , 
h “ hjj - vt 


and the nmv Hamiltonian la 


F - . vll . 


UUU*t.s.4- 


UIK 

lu I'iHJR 


The equations «tf motion are then 


dL 

8F 
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8F 
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3F 
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dt 
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with F - Fq f Fj * Fg: 
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li 


0 ■ 2 ' ’ 
2L 


I'l -•’21*7 




1^2" E\m ■ 


in the above, u - g t f , and f is the true anomaly. 


(12) 


(13) 
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The Ilamiltwiion la now expanded In tema of powers of tlio small parameter J^; 

I'y OtJg) , 

Fg - 0(J^) . 

2. 2 First-Order Perturbations 

To simplify the writing, we shall now consider only one value of kj, as defined 
In Section 2. 1, jmd drop tlie index in expressions such as DCk^) and The 

Hamiltonian is then such that 

Fg - D cos 20j^ . 

We define 0 by 



To derive the complete perturbations, wo shall sum the individual expressions obtained 
for each value of kj. 

In the process of removing short-period terms from the Hamiltonian in I’aper I, 
wo applied Hori's method by lie series (Ilori, 19GG) and considered a canonical 
transfonnation 


(L,G,H,^g,h) 

F 


(L',GMr,i",g',li') 

F' 


such that 


I. - 1/ + ^ t- 


i .1- qI , 

r • r' I * i JM. , 

9g' 2 U)g'’®i "• ’ 


} 
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II = H' * 


i - i' - 
K == g' ' 

li V. h' - 


as 


W 



8G' 



ns 1 |.')S 

m' ”2 


> 


} 


} 


} 


as <lescrll)6d by equatloiis (30) in Paper I, 


The transformation Involves a generating fuisction 




(H) 


where 


SR = Z, ! 7u 


1 


^1/2 perturbations due only to resonant harmonics, SJ the perturba- 

tions due only to Jg (or to the aonals, if otlier zonals arc also considered), and Zg 
th/fi; perti’.rbations due to the interaction of Jg with the resonant tesserals. 


bi order to avoid small divisors, S is developed in terms of powers of the square 
root of the small parameter Jg, as is tlie custom in resonance problems (Paper I). 

We reintroduce the characteristic small parameter of the problem: 


Y= V 



To find Sj^g, 
(Hori, 1966) 


Z^, and Zg, we gather terms of the same order relative to Jg from 


F'^ F + {p,s} + |{{r,s},s}+ ... , 


(15) 


whei-e F' is the now Hamiltonian, and y is assumed to be of order 1/2. 
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In comparing this theory with numerical integration, we found good agreement for 
tlie perturbations of highest order due only to the resonant tesserals. The relevant 
part of the determining function is S Appendix A, we repeat the derivation of 
expressions of the derivatives of this function with respect to the modified 
Delaunay variables, and definitions of the quantities A, Ig, Iq, k, and A. 

Some problems developed, however, when the Jg tesseral-interactlon terms 
were considered. It turns out (cf. Paper I) tlrnt the determining function 7,^ cannot 
be defined by 


yA 




+ b 


^ 1/2 
2 9g' 


+ b 


^1/2 ) 

3 W / 


) 


where 


b 


1 



8Fi 

^2"" W 





} 


since here we have neglected the contribution of tlie secular rate to the variable g. 

To the highest order, then, let us consider instead the following equation as a 
definition of Zgt 

{FQ + Fj,,Zg} + {F^,Sj^/g}=0 . (16) 

We assume that Zg is a function only of L', G', H', and G, where we recall 

20 = al' + (3(g' + h') + ip(g') - IT . 

Then we obtain 
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azg 





9^2 

80 » 
W “ 2 80 » 


wltii 


cl = p + 


ii 

8g' 


f 


which yields 


gy f g ^^2 
2 80 


^^l’®l/2^ = •"O'CgY - ^2^4°2^ ^^2 ’ 


where 


<J - bj^g + bgd + bgP 


and 


( 17 ) 


(18) 


n 1 ® 

^4 “ D 8g 


[cf. equation /44) of Paper I]. Hence, 

Z2=Ba2~Q) , (19) 

where 

2gc„y 

B= ^ . 

gy + a 
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The derivatives of Zg with respect to the modified Delaunay variables are as 
follows. If X is a modified Deli.unay variable, then 


ax ax ^^2 


- 0) f B 


M" ■" ax/ » 


( 20 ) 


where aig/ax and 30/8X are given in Appendix A, and 


^ A g \ . c 
ax" "av+g ax \; "aY + gy 


2g 

a\ + g 




a 


ay + g/ 


J 


in which 


e= 1 

if 

X = L' 

!l 

o 

if 

X L' 


and 


r._ 8b, Sbf, abo ^2, 

.221 n — i + d — - 4. 3 2 4. h .2-2 

ax ax “ ax ax ”2 ag ax ’ 

The relevant expressions are as in equations (44) and (65b) and Appendix A of Paper I; 
they are repeated in Appendix B herein. 

2. 3 Second-Order Peiturbations 


Instead of considering the second-order perturbations due to the resonant 
tesserals separately from those due to the interaction terms, let us define a determin- 
ing function Zg by 

{Po + Fj,Z3}+|{{Fj+Fj,Sj/2 + Z2}, Sj^2-22}=» ■ PD 

This yields, assxuning Zg = Zg(L', G', H', 0), 
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(22) 


°"^2 80 ^ ^P®l/2 ^1/2 ^ 2 } » 


which, after some algebra, gives 


8Z / \ 2 2 

3 « / . 1 \ ^ „ r X k sin 0 cos 0 , „ „ k sin 0 cos 0 

Rj, (A 1 (A-1) -f Rg ag-V A ^3 <^2-®> A 


+ R, (A-1)" , 


(23) 


where 


R^ = R2=-(aC2Y“) 


2.2 


X 


^3-- 


-aC2Y^(“^ + 


4aCj^Y 


+ 1 


(24) 


R 


2„ .3 

4 ’ aY + cr 


(3aY + 2a) 


In equations (24), terms of higher order were neglected and the following were used; 


X- 


aa^^ + dag + pUg - PgU^ 


0 =-ii 

.2 


3|j a 


c„ = 


2 „ 2 ’ 
3|i a 


o =i-^ 

D 0L' ’ 


= 1 ® 

^2 D 8G' ’ 


= 1^ 

^3 D 8H' ’ 


•n _ 3B , J 8B , n 8B 
®p 8L' ^ 8G' 8H' • 




OB'* THE 

, VOiM 
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After some algebra, then, we obtain 


r 


2^3 = ^! 


1 - 0 + ^ sin 20 


- Pgl2 + P3O + I’4lo > 


( 25 ) 


where 


= 2Rj^ + R^ + Rg , 

Pg = R^ + 2R^ + Rg + ^(Rj + Rg) , 
Pg = R^ - 2R^ + ARg , 

P4=Hj(1 + A) . 


The derivatives of Z„ with respect to the modified Delaunay variables are as follows; 

o 


aZg 0Pj 

W 


K 


1 2\ 1,2 

1 1A 1 *> OA 

1 ~ S JU -io 


ax ‘2 ' ax ° ax "O 


0 I 9 on 91f) 

^ p 2 + P — + P — + P ,0 

^2 dX 0X 4 8X ' 


(26) 


where 


00 






2 2 


2 “x 


sin 20 


In equation (26), the derivatives aig/aX, 00/0X, and 8 Iq/ 0X are given in Appendix A. 
Now we get 


0R 


i = - (aCgY^)^ ^ (4ttC2Y^) (1 + 4aCj^V) > 


dRr 

c 

w 


= £ 


~2y (2ac^y + 1) + 4aC + ly - 4aCgY^ 
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®4 

a3T 


3 4 

aCgY 


(ay + O') 


8a 

IW 


+ 6 


(- + 8c ,YCr + ^ i- 


oy 


a ay + a 


where 


£>= 1 if X- L' , 

£ « 0 if X ^ L' . 


In conclusion, the determining function for tlio canonical transformation used to 
remove short-period terms from the Hamiltonian is 

S = Sj^/2 + Zg + Zg , 

and the new Hamiltonian is 

F'= Fjj + F'^ + F' , ( 27 ) 

where 

Fj(L',H')=Fo(LyHO , 

F'(L',G',H')-Fj(I/,G',H') , 

F' (L', G', H', g') = - X) • 

The new equations of motion are then 


dL' 

dt' 

8£' ^ » 

II 

ife 

8F' 

~ W 

dG' 

_ 8F' 

dg'_ 

8F' 

dt' 

3g' ’ 

dt' “ 

" 8G' 

dH' 

9F' . 

dh' 

0F' 

W 

8h' " ° ' 

dt' " 

■ W 


and the long-period variable g' remains to be removed, following the procedure given 
in Section 5 of Paper I. 
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3. COMPARISON OF THE THEORY WITH A NUMERICAL INTEGRATION 

The theoretical computations have been compared with a numerical Integration 
provided by the routine DVDQ. Tlie satellite considered was Geos 1, which passes 
through the a = 1, p = 12 resonance with the following orbital elements; 

semimajor axia Uq = 8. 074 X 10 m , 
eccentricity e^ = 0. 073 , 

Inclination I^- 59® . 

We take tlie constants of the problem to be 

^0 ” ’ 

^0" ^0 ’ 

IIJ^GgcosI , 

(29) 

f5=0 , 
gg=0 , 
h"=0 . 

2 -4 

In tills case, k is of the order of 10 , so the motion of tlie satellite is of the circula- 
tion type, The period of the short-term perturbations is approximately 7 days. 

After removing the short- and long-period perturbations, we obtain errors in 
the Delaunay variables L, G, and H for a period of 8 days, witli computations per- 
formed every 0. 5 day. The results are shown in Table 1. 
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Table 1. Errors I’"* uction variables as computed over a period 
of 8 dsy j. 


Variable 

Amplitude of ^ 
perturbation X 2 

Error 

Relative error 
(%) 

L 

48 X 10“^ 

9. 75 X 10"’^ 

0.2 

G 

67X 10"® 

8. 86 X 10"® 

13 

H 

57 X 10"® 

1. 14 X 10"® 

0.2 


^ A 

The semimajor axis is ejqpressed in units of 10° m. 


We obtain the following errors in the Delaunay angular variables: 

_7 

i - 7. 79 X 10 rad in 8 days, corresponding to '^0. 8 m per day in the 
ascending node of the orbit (a= 8. 074 X 106 m), 

_7 

g = 7. 05 X 10 rad in 8 days, corresponding to 0. 75 m per day in the 
argument of perigee (u = 8, 074 X 10® m), 

h = 8. 7 X 10"^ rad in 8 days, corresponding to 0. 002 m per day in the 
satellite's position in the orbit (a= 8. 074X 10® m). 

Moreover, the' errors in f and g have opposite signs, so that the error in f + g is of 
—8 —8 

the order of 8 X 10 rad in 8 days, or 10~ rad per day (see Figure J). 

We conclude that there is satisfactory agreement in all variables except G. 
Adding or removing the contributions of second-order perturbations (the generating 
function Zg) does not affect tlie results. We suggest a tentative explanation for the 
large error in G. In the process of removing the variables St and h from the 
Hamiltonian, we have written 


where 
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We took Zg from equation (4), whore we asBxuned that Zg depends on g' only through 
the rescaiant variable 0, However, after obtaining Zg and its derivatives, we get 
[see equation (20)] 


0Z 


2_ 8B 


^ “ 8g' '^2 


(lo-O) f-B 




1 ) I 


a 


2 ^ 2 ' 


'V 


(30) 


where d~ (3 + 4'(s') 8-4 ■' (fA^) (8k^/dg') are functions of L', G', H', and g', and 

both have significant derivatives witli respect to g'. Hence, there is an Inaccuracy in 
the dependence of Zg on g'. The correct equation for Zg is 



8F, 0Z„ 

■g^ “ --{FpSi/g} , 


(31) 


where 


Zg= Zg(I.,G,H,0,g) 
and 

dg'’" “ 00 0g' 8g^ 

Another factor in support of this explanation is that, when Jg is set equal to zero 
(and, hence, OFj^/OG' ”■ 0), the error in the variable G is of the order of 0. 2%, as it 
is for the variables t and II. 

It seems that, to obtain better agreement in the variable G, we must take into 
account the intricate dependence of Zg on g' and solve equation (31) more accurately. 
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APPENDIX A 


EXI’HESSIONS OF S j/2 ^ND PI’S DMIIVATIVES ^ITII KESPECT TO THE 
MODIPIED DEIAIINAY VAHIABLES 


The ftmctlon Sjy2 obtained by Holving the reBcmance equation 


where 


il v-ii— 

,3 


Y - f V 


and 


If 


«S,/y 


A- gyJ'''* 

^ 2 - > 


3}i akj^ 


and 


,2 12 1) 

T 2 > 

1/ Y 


A'^ Vl - sin^ 0 , 


(A-1) 


A-1 



then, solving equation (A-1) and applying some ccntinalty arguments yield 


33 

W 


1/2 _ 


A(A- 1) 


Ri'T‘i:'v ; f;v or tiik 

Oklijii.'.AL J,’IM/li 


and 


^ 1/2 ” '^^2 "" ’ 
where we define 

l2=E(G,k) . 

Let us also define 

Io=F(0,k) . 

Here F(0,k) and E(0, k) are the elliptic integrals of the first and second kind, respec- 
tively. The derivatives 0^82/2 

®^l/2 

a^-=aC2Y(A-l) , 

®®l/2 

■gljT ~ dCgY (A-1) + > 

®^l/2 

(A-2) 

®®l/2 2 

gx7^= f (Iq-®) SCfY (Ig-O) + C2Pl'^(^"^^ ^^2"V ’ 

go? = Cg Y P2 (^ " 1) + °2 ^ ^2 ~ ^0^ ’ 


^® l/2 

3H' ~ V Cg P3 (A - 1) + Cg Y ^3 (I2 “ ^0^ ’ 


A-2 





where Cp Cg, ag, are as given in equations (24), 


1 0D 

^ 4 "J 5 W ^ 


P 


_ Bi!> 

1 0L' ’ 


^2 0G' ’ 


d4> 


3 0II 


/ ) 


■4 0g 


Wo also have 


80 _Pl 


80 _ a 

W 2 » 


00 _ ^2 
3G' " T > 


80 




(A-3) 


_3 

8II' ■“ 2 » 


30 

W 


and, if X is a modified Delaunay variable, 


0A_ Sc f ^ 1 \ sin 0 cosO 00 

OX 2 V^"Ar A W > 


Of 


a., 


^ ^ A 1 . ^ /T T \ 

3X OX ^ 2 '^2 " ■^O^ » 


(A-4) 


where 


^ = 190. + !x/t 

3X A 3X 2 ^-2 ^0^ > 


u 

W0>"j 


1 -k 


dx „ 1 f T K sin 0 cos 0 

2 VZ" A 
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In the above, 



if X^-i'orh' 
if X= 1/ 

if X= G' 

if x==ir 
if X = g' . 
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APPENDIX B 


DERIVATIVES OF bg, bg, AND ■4>(g'), AND SECOND 
DERIVATIVES OF l'(g') 

We recall the definitions of D and 0, introducing the quantities B and C: 
—2 ^2 

D ” B + C , 

D cos ip - B , 

D sin C . 

Then, if X is any of the variables L', G', IF, or g', we get 


^ fOX “ ^ ax ^ 8X ' 


(B-1) 


If X g', we have 


kiP-i-ko 


m 

ax' 




kg x=0 




dC 

ax' 


k^(3+kg 

z i: 

kg x=0 


8_ 

ax 


S(kp kg) A(kp kg, x)J sin [(kg - 2x)g - 


(B-2) 


°^ yi' ''O’ !!b. ”upj^ ^ Ofe) X, °°n,p,q<^> ae 

ax ai ax^'®^ '^n,m,p'^^ ae «'5^ ’ 


ax 


as(kjkg) _ ) 0 

ax "j 2 (pkj+kg + i) 


if X ^ L 


S(kpkg) if X= L ; 
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and if X =• g', we get 


Pkj t Uq 

W \ 

1<Q x=0 ^ 

pkj fkQ 

li XI 2 ®^^^P ^^0^ ^<^^P [<^^0 " " ^k 

ko ^‘=0 

I*'or the derivatives of i{i, we have 
C = D sin ‘4‘ ) 


so that 


dC 

ax 


= D 



cos !/» + 


ax 


sin ii 


and 


^_i./ac l^p\ 
ax “ ' B \ax " D ax 7 • 


We also have 


ab, abn abn 

“ X=i'orh'org' , 

8b, 4b, 

"17 » 

ab2_ Sbg 

aI7 " XT' > 

■^ " L' ' 

8b^ 3b ^ g Jgl/ll^n'^ 

_ = _ 4 . _ ______ ^ 
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''J ’ 

(B-3) 

'Kq] ■ 


(B-4a) 

(B-4b) 


(B-5) 



Obg 4b 

A. 9 ^ 

i^jgi/ii-^ii- 

2 jy3j.,7 

abg f)bg 
W" ‘"W 

} 

0 J 

2ii‘‘n^H' 

mr" "2 

L/4 q/H ^ 
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0bg 3 
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